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1 The probe method for inverse obstacle scattering 
problems at a fixed wave number 

In this paper we consider inverse problems for partial differential equations. We restrict 
ourself to the reconstruction issue of the problems and refer the reader to [29] for several 
aspects and uniqueness results in inverse problems for partial differential equations. 

More than ten years ago Ikehata discovered two methods for the purpose of extracting 
information about the location and shape of unknown discontinuity embedded in a known 
background medium from observation data. The methods are called the probe and enclo- 
sure methods. This paper presents their past and recent applications to inverse obstacle 
scattering problems of acoustic wave. 

The probe method was originally introduced in 1997 and published in [5] . Since then 
the method has been applied to several inverse problems for partial differential equations 
[6, 7, 13, 19, 20] and still now some new knowledge on the method itself added in [16, 21]. 

In this section we present one of typical applications of the probe method published 
in [7]. Therein the author considered an inverse obstacle scattering problem at a fixed 
wave number. We denote by D and B R an unknown obstacle in R 3 and open ball radius 
R, respectively. We assume that: D is an open set with smooth boundary satisfying 
D C B R and that B R \ D is connected. dB R indicates the location of the emitters and 
the receivers. 

Let k > 0. Given y G dB R let <&(#) = $D(x,y; k),x G R 3 \ D denote the solution of 
the problem: 

(A + A; 2 )$ + 5( ■ - y) = in R 3 \ D, ^ = 0ondD 
and the outgoing Sommerfeld radiation condition lim r(d§/dv — ik<&) = 0, where r = \x\ 

r — >oo 

and v is the outward normal relative to D. 

Inverse Problem 1.1. Fix k. Reconstruct D from the surface data $d(£, y; k) given at 
all x G dB R and y G dB R . 

The $d has the form &d(x, y, k) = $o(^, V, k) + E D (x, y; k), where E(x) = E D (x, y; k) 
satisfies 

(A + A; 2 )E = 0inR 3 \A ^ = ondD 

ov ov 

and the outgoing Sommerfeld radiation condition lim r(dE/du — ikE) = 0; &o(x, y; k) = 

e %k \ x ~ y \ / (4tt\x — y\). The E D (x,y;k) is called the scattered wave field generated by the 
point source 5( - — y) located at y. $£>(x,y; k) is called the total wave field. 
In [7] the author has established the following result. 

Theorem 1.1. Assume that k 2 is not a Dirichlet eigenvalue for —A on B R nor an 

eigenvalue for —A on B R \ D with homogeneous Dirichlet boundary condition on dB R 
and Neumann boundary condition on dD. Then one can reconstruct D from V\ k) 

given at all x G dB R and y G dB R . 

A brief outline of the proof is as follows. Set Vt = B R . We starts with introducing two 
Dirichlet-to-Neumann maps for the Helmholtz equation in f2 \ D and Q. 

Given / G if 1/,2 (<9f2) let u G H 1 (Vl\D) be the weak solution of the elliptic problem 

(A + A; 2 )m = ian\D, — = 0ondD, u = fondtt. (1.1) 
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The map Ad : / i — > du/dv\gn is called the Dirichlet-to-Neumann map associated with 
the elliptic problem. Set also A D = A for D = 0. 
Theorem 1.1 is divided into two steps. 

Step 1. One can calculate A — Ad from $d(x, V\ k) given at all x G dQ and y G dQ. 
Step 2. One can reconstruct D itself from the integral / (A — A D )f • fdS for infinitely 

J dQ 

many fs independent of D. 

Note that the integral in Step 2 has the form 

/ (A -A D )f-JdS=[ (^u-^v)dS 
Joci Jan \dv dv J 

where v = v{x), x G Q solves (A + k 2 )v = inf2, v = f ondQ; u = u(x), x G Q\ D 
solves (1.1) with / = v\qq. Thus infinitely many f means infinitely many v. 
The step 1 consists of two parts. 

(i) Given / find the solutions g and h of the integral equations 

/ $ (x, y; k)g(y)dS(y) = f(x), f $ D (x, y; k)h{y)dS{y) = f(x), x G dQ. 
Jan Jan 

(ii) Compute (A — A D )f by using solutions g and h in (i) by the formula (A — A D )f = 
g-h. 

This type of procedure, like (i) and (ii) have been known for the stationary Schrodinger 
equation [35] and the proof is an adaptation of the argument. Thus the point is Step 2. 

1.1 Step 2. 

In this subsection we explain Step 2. Instead of the original formulation of the probe 
method we employ a new one developed in [16, 21]. 

1.1.1 Needle, Needle sequence 

Definition 1.1. Given a point x G Q we say that a non self-intersecting piecewise linear 
curve o in Q is a needle with tip at x if a connects a point on dQ with x and other points 
of a are contained in Q. We denote by N x the set of all needles with tip at x. 

Let b be a nonzero vector in R 3 . Given x G R 3 , p > and 9 g]0, 7r[ set C x (b, 9/2) = 
{yE R 3 | (y - x) ■ b > \y — x\\b\ cos(#/2)} and B p {x) = {y G R 3 \\y - x\ < p}. A set 
having the form V = B p (x) H C x (b, 9/2) for some p, b, 9 and x is called a finite cone with 
vertex at x. 

Let G(y) be a solution of the Helmholtz equation in R 3 \ {0} such that, for any finite 
cone V with vertex at 

/ \VG{y)\ 2 dy = oo. 

Hereafter we fix this G. 

Definition 1.2. Let a G N x . We call the sequence {v n } of iJ 1 (fi) solutions of the 
Helmholtz equation a needle sequence for (x, a) if it satisfies, for any compact set K of 
R 3 with K C Q \ a 

JimJIM • ) - G( ■ - x) \\ LHK) + || VK( •)-<?(•- x)}\\ LHK) ) = 0. 
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The existence of the needle sequence is a consequence of the Runge approximation 
property (cf. [30]) for the Helmholtz equation under the assumption on k: k 2 is not a 
Dirichlet eigenvalue for — A on Q. See the appendix of [7] and A. 1. Remark in the appendix 
of [16] for the proof. The unique continuation property of the solution of the Helmholtz 
equation is essential. 

1.1.2 Special behaviour of the needle sequence 

In the following we do not assume that k 2 is not an eigenvalue for — A in Q with Dirichlet 
boundary condition. 

Lemma 1.1. Let x E Q be an arbitrary point and a G N x . Let {v n } be an arbitrary needle 
sequence for (x, a). Then, for any finite cone V with vertex atx we have \\ Vv n \\i,2(y^ — > 
oo as n — > oo. 

Lemma 1.2. Let x E Q be an arbitrary point and a G N x . Let {v n } be an arbitrary 
needle sequence for (x,a). Then for any point z G a and open ball B centered at z we 
have || Vi>„||i,2( Bn Q) — y oo as n — > oo. 

Note that from Definition 1.2 and Lemmas 1.1 and 1.2 one can recover a G N x itself from 
the behaviour of any needle sequence for (x, a). 

Summing up, we see that {v n } has two different sides: 

(A) converges to singular solution G(y — x) with singularity at y = x outside a; 

(B) blows up on a. 

These different sides of needle sequences yield two sides of the probe method which we 
call Side A and Side B. 

1.1.3 Indicator function and Side A of the probe method 

Let v satisfy (A + k 2 )v = in Q and u solve (1.1) with / = v\qq. Set w = u — v inQ\D. 
The w satisfies 

(A + k 2 )w = 0inn\D, w = ondtt, — = -— ondD. (1.2) 

ov ov 

Integration by parts yields 

/ (A $ -A D )(v\ dn )-vdS= [ \Vv\ 2 dy-k 2 f \v\ 2 dy 
Jan Jd Jd 



+ [ _\Vw\ 2 dy -k 2 / _\w\ 2 dy. 
Jn\D Jn\D 

This motivates 

Definition 1.3. The indicator function I(x), x G Q \ D is defined by the formula 

I(x)=[ \VG(y - x)\ 2 dy - k 2 f \G(y - x)\ 2 dy + / _\Vw x \ 2 dy - k 2 f _\w x \ 2 dy, 
Jd Jd Jn\D Jn\D 

where w x is the unique weak solution of the problem: 

(A + k 2 )w = infi \ D, — = -tt(G( ■ ~ x)) onD, w = ondQ. 



1.3) 
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The function w x is called the reflected solution by D. 

The following theorem is based on the convergence property of needle sequences and 
says that 

• one can calculate the value of the indicator function at an arbitrary point outside D 
from A — A D ; 

• the indicator function can not be continued across dD as a bounded function in the 
whole domain. 

Thus one can reconstruct dD as the singularity of the field I(x) which can be computed 
from the data with needles and needle sequences. That is the meaning of the following 
result. 

Theorem A. It holds that 

• (A.l) given x G f2 \ D and needle a with tip at x if a D D — 0, then for any needle 
sequence \v n \ for (x,a) we have I(x) = lim / (An — ^D)(vn\dn) -v^dS; 

• (A. 2) for each e > sup {I(x) \ dist(x, D) > e} < oo/ 

• (A. 3) for any point a G dD lim I(x) = oo. 

The key for (A. 3) is to establish lim sup ll^ll^^m < oo. An outline of the proof is 

x — >a 

as follows. Using the solution of the boundary value problem: (A + k 2 )p = w x in f2 \ D, 
p = on dfl and dpjdv = on dD, we have the expression 



J \wx\ 2 dy = J 9d (p( x ) ~ P(y))~^-(y - x)dS(y) + k 2 p(x) $ (y - x)dy. 

Applying a standard regularity estimate of p: \\p\\h^{q\d) — C \\ w x II L 2 (n\D) anc ^ ^ ne Sobolev 
imbedding: \p(x) - p(y)\ < C\x - y\ 1/2 \\p\\ H 2 {n \D), x,y G Q \ D and |b|| L oo (n \S) ^ 
C|bll/f 2 (Q\D) ^° ^ s right-hand side, one gets an upper bound of H^H^^yD) which in- 
volves integrals of weakly singular kernels over dD and D. 



1.2 Remark I. Side B of the probe method and an open problem 

Since mathematically Theorem A is enough for establishing a reconstruction formula, in 
the previous applications of the probe method we did not consider the following natural 
question. 

• Let x G VL and a G N x . Let £ = {v n } be a needle sequence for (x, a). What happens on 
the sequence 

I(x,a,£) n = / (Ao - A£>)(u n |an) -v^dS,n = 1,2, • • • 
Jan 

when x is just located on the boundary of obstacles, inside or passing through the obsta- 
cles? We call sequence {I(x, a, £)„} the indicator sequence for (x,a) and £. 

In practice the tip of the needle can not move forward with infinitely small step and 
therefore in the scanning process with needle there is a possibility of skipping the unknown 
boundary of obstacles, entering inside or passing through obstacles. So for the practical 
use of the probe method we have to clarify the behaviour of the indicator sequence in 
those cases. The answer to this question is 
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Theorem B. Assume that k 2 is sufficiently small (not specify here). Let x G f2 and 
a G N x . If x G Vi\D and aHD ^ or x G D, then for any needle sequence £ = {v n } for 
(x, a) we have lim I(x, a, £) n — oo. 

In the proof the blowing up property of needle sequences is essential. 

A sketch of the proof For simplicity, we consider here only a single obstacle case. We 
make use of two well known Poincare's inequalities: 

(I) < C(n\D)\\Vw\\ 2 for all w G \ Z>) with w = on dft; 

(II) \\v - v D \\ 2 L 2 {D) < C(D)\\Vv\\ 2 L2[D) for all v G if^D), where w D = J vdy/\D\. 

Let A be an arbitrary Lebesgue measurable set with A C D, \A\ > and t> G L 2 (D). 
A simple argument in [42] gives ||i> — fA||L 2 (D) — 2iiCi||i; — itoll^p), where va — J vdy/\A\ 
and K A = 1 + A combination of this and (II) yields 

j D \v\ 2 dy < AK A C(D) j D \Vv\ 2 dy + 2\D\\v A \ 2 . (1.4) 

Let u = u n solve (1.1) with / = v n \sn and set w n = u n — v n . It follows from (1.3), (I) 
and (1.4) that 

I(x, a, On > (1 - k 2 C(Q \D))[_ \Vw n \ 2 dy 

JQ\D 

+ (l-4k 2 K A C(D)) I \Vv n \ 2 dy-2k 2 \D\\(v n ) A \ 2 . 

JD 

Thus if k satisfies k 2 C{Vt \ D) < 1, then we have 

I(x,or,0n>(l-4k 2 K A C(D)) ( \Vv n \ 2 dy-2k 2 \D\\(v n ) A \ 2 . 

JD 

Write 1 - Ak 2 K A C{D) = 1 - 8k 2 C(D) - Ak 2 (K A - 2)C(D). Here we make k smaller in 
such a way that 8k 2 C(D) < 1. Using an exhaustion of fi \ <r, one can construct A G D 
in such a way that |A| w |D| and A C \ u. Since — 2 = — 1, one gets 

1 — Ak 2 K A C(D) > 0. Note also that the sequence {{v n ) A } is always convergent for a fixed 
A. Thus the blowing up property of the indicator sequence is governed by that of the 
sequence {|| Vv n \\ 2 L ^ D) }. 

A combination of Theorems A and B yields another characterization of the obstacle. 

Corollary 1.1. Assume the smallness of k 2 same as Theorem B. A point x G Vt belongs 
to Q\D if and only if there exists a needle a with tip at x and needle sequence £ for (x, a) 
such that the indicator sequence is bounded from above. 

Needless to say, this automatically gives a uniqueness theorem, too. 

An open problem in the foundation of the probe method is the following. 
Open problem 1.1. Can one remove the smallness of k 2 in Theorem B? 

Here are some closely related technical questions. 
• Is it true ?: if x G f2 \ D and a C\ D ^ % or x E D, then 

lim „ IKII f (D) = 0. (1.5) 
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• Let u = u n solve (1.1) with / = v n \gn and set w n = u n — v n . We know that if x G D, 
then || Vu> n || L 2(^\75) — > oo as n — > oo ([20]). The question is: identify the points in 
Q\D that really contribute the blowing up of Vw n . See [16] for an example in the case 
when k — 0. 

• Is it true ?: if x G VL \ D and (inD^OoriGD, then 

hm = 0. 

n^oo ||Vw„||l2(D) 

See [16, 19, 20] for more information on these questions. 



1.3 Remark II. An explicit needle sequence 

From Lemmas 1.1 and 1.2 we know that given a G N x the energy of an arbitrary needle 
sequence {v n } for (x, a) blows up on a. However, it will be difficult to understand the 
behaviour of v n (y) at each y G a. In this subsection, we give a family of special solutions 
of the Helmholtz equation with two parameters that yields an explicit needle sequence for 
a straight needle. We call such a family a generator of needle sequence. 

The contents of this subsection are based on the classical materials developed by 
Yarmukhamedov, Mittag-Leffler and Vekua. 



1.3.1 Yarmukhamedov 

The following fact is taken from the article [45] . 

Theorem 1.2. Let K(w) be an entire function such that: K(w) is real for real w; 
K(0) = 1; for each R > and m = 0, 1, 2 sup (w)\ < oo. 

\Rew\<R 

Define 

c 2^ , s f°° t ( K(w)\ du 
Jo 



w 



\J\x'\ 2 + u 2 

where w = X3 + ^l^'l 2 + u 2 and x' = (xi,x 2 ) 7^ (0,0). Then one has the expression 
&k(x) = l/(4vr|a:|) + Hk(x) where Hk satisfies AHk(x) = mR 3 . 

Note that <3>^ can be identified with a unique distribution in the whole space and 
satisfies A$ K (x) + 5(x) = inR 3 . 

Example 1. K(w) = I. In this case we have $^(x) = l/(47r|x|). This is because of 



-c 



du 



Ak\x\ J-00 4tt 2 (\x\ 2 + u 2 

and 

1 

— Im 



\X\ 2 + U 2 y x3 + i ^\ x >\2 + u 2 J ^ X l\2 + U 2 

Thus for general K we have 

1 r°° (K{w)-l\ du 
Hk(x) = -— 2 J o Im 



w 



sj\x'\ 2 + u 2 
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Example 2. K(w) = e TW . r > a parameter. In [12] the author pointed out that <&k(x) 
with this K coincides with the Faddeev Green function G z (x) with z = r(e 3 +iei): 

G *( x ) = tttts / — ^ — dr i- 

(27r) 3 Jr3 |r/| 2 - i2z ■ i] 

The Faddeev Green function has been applied to several inverse boundary value/scattering 
problems by Sylvester-Uhlmann [43], Novikov [36], Nachman [35], et al.. 



1.3.2 Mittag-Leffler 

Let < a < 1. The entire function of the complex variable w 

w w 2 w 3 

W = + r(l + a) + r(l + 2a) + r(l + 3a) + " ' ' 

is introduced in [34] and called the Mittag-Leffler function. 

It is known that K(w) = E a (rw) with r > satisfies the condition in Theorem 1.2 (cf. 
[2]). In [46] Yarmukhamedov applied this function with a fixed a to the Cauchy problem 
for the Laplace equation in two dimensions. 



1.3.3 Vekua 

The Vekua transform v \ — > T k v in three dimensions [44] takes the form 

T k v(y) = v(y) - ^ v(ty)Ji(k\y\y/T^t) sf^dt 

where J\ stands for the Bessel function of order 1 of the first kind. 

The important property of this transform is: if v is harmonic in the whole space, then 
T k v is a solution of the Helmholtz equation Au + k 2 u = in the whole space. 



1.3.4 Generator of needle sequence 

Using materials introduced by Yarmukhamedov, Mittag-Leffler and Vekua, the author 
found an explicit needle sequence when the needle is given by a segment. 

Given < a < 1 and r > define v(y;a,r) = —H K (y), y G R 3 , where K(w) = 
E a (rw). This v is harmonic in the whole space and thus the function v h (y;a,T) = 
T k v(y; a, r), y G R 3 satisfies the Helmholtz equation in the whole space. 

Theorem 1.3([21]). Let x G f2 and a be a straight needle with tip at x directed to 
u = (0,0, 1) T , that means: o has the expression o = {x + s u [0 < s < 1} with I > 0. 
Then the function v k ( ■ — x; a,r)\n as a — > and r — > oo generates a needle sequence 
for (x, a) with G = G k given by 

( Jk\y\ \ 



Note that since the function 



sin k\ 


y\ 


4ir\ 


y 





,yeR 3 
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satisfies the Helmholtz equation in the whole space, the function 



sin k 


y-x\ 


4tt; 


y-x\ 





v [y — x; a,r) + % — — ; — , y G \l 



generates also a needle sequence for (x, a) with 

p ik\y\ 

G(y) = -— v (1.6) 
47r|y| 

Thus now we have an explicit generator of a needle sequence for a straight needle with 
(1.6). This makes the probe method completely explicit in the case when one uses only 
such a needle. Everything is reduced to the choice of small a and large r. 

This is very important also in the singular sources method by Potthast [38] since in 
his method one has to construct the density of the Herglotz wave function (cf. [3]) that 
approximates locally the fundamental solution of the Helmholtz equation in a domain like 
Q\o~. However, Theorem 1.3 shows that instead one can consider only a simpler problem: 
construct the density of the Herglotz wave function that approximates v k (y — x; a, r) on 
the whole boundary of a geometrically simpler domain like a ball. 

Open problem 1.2. It would be interesting: do the numerical testing of the probe and 
singular sources methods in three dimensions with this explicit needle sequence. 
Open problem 1.3. A mathematically interesting question is: find a generator of a 
needle sequence for a general needle. 

Note that Yarmukhamedov [47] made use of ^x(y — x) itself not its regular part 
Hxiy — x) to give a Carleman function which yields a representation of the solution of 
the Cauchy problem for the Laplace equation in three dimensions. 

Finally we give a remark that is closely related to Open problem 1.1. In [21] an explicit 
formula of the precise values of v k (y — x; a, r) on the line y = x + su(—oo < s < oo) is 
given. They are: 

• if y — x + sou with s ^ 0, then 

v\y - x- a, r) = ^ ^M~ cos ks _ A f\i _ w *)- 1/2 E a (T(l - w 2 )s)J l (ksw)dw; 

4:71 S 47T JO 

• if y — x, then 

v k (y - x;a,r)|j 



ly ~ x 47rr(l + a)' 

Moreover, we see that Vv k (y — x; a, r) on the line y = x + su (— oo < s < oo) is parallel 
to oj. In particular, we have 



r 2 



v, (v-*;<*,T)\ v=x = Anni + 2Q) u. 

It seems that the behaviuor of v k (y — x; a, r) and its gradient at y = x suggest the validity 
of (1.5). 
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2 The enclosure method for inverse obstacle scatter- 
ing problems at a fixed wave number 

The enclosure method was introduced by the author in [10] and has been applied to 
several inverse problems for partial differential equations. In this section we present its 
applications to inverse obstacle scattering problems at a fixed wave number. 



2.1 The enclosure method with infinitely many data 

The method applied to inverse obstacle scattering problems is based on the asymptotic 
behaviour of the function (we call the indicator function again) 



t\ — yf (A - A D )(v\gn) -vdS, 



where v = e x < TU1 + i ^ T2 + k2ul± ) having large parameter r; both oj and oj 1 are unit vectors 
and perpendicular to each other. 

This v satisfies the Helmholtz equation Av + k 2 v = in the whole space and divides 
the whole space into two parts: if x • oj > t, then e~ rt \v | — > oo as r — > oo; if x • oj < t, 
then e~ rt \v | — > as r — > oo. 

The method yielded the convex hull of unknown sound-soft obstacles by checking the 
behaviour of the indicator function. It virtually checks whether given t the half space 
x ■ oj > t touches unknown obstacles. 

In [9] an extraction formula of an sound-hard obstacle D C R 3 with a constrained on 
the Gaussian curvature of dD from Dirichlet-to- Neumann map A^ has been established. 
Its precise statement rewritten with the present style is the following. 

Let us recall the support function of D: hoioj) = sup x ■ u, u G S 2 . The convex hull 

of D is given by the set ^ & s 2 { x G R 3 | a; • w < hc(u)}. Therefore, knowing h D {oj) for a 
oj yields an estimation of the convex hull of D from above. 

Theorem 2.1. Assume that the set {x e dD \ x ■ oj = h£,{oj)} consists of only one point 
and the Gaussian curvature of dD doesn't vanish at the point. Then the formula 



lim — log 

is valid. Moreover, we have: 
if t > hoioj), then 



(A - A D )(v\ dn ) ■ vdS 
an 



h D (u), 



lim / (A - A D )(e- Tt v\ dn ) ■ e-^vdS = 0; 

if t < hr>{uj), then 

lim / (A - A D )(e~ Tt v\ dn ) ■ e~ rt vdS = oo; 

if t = h D {oj), then 

liminf / (A - A D ){e- Tt v\ dn ) ■ e^dS > 0. 
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Note that: if one considers the Dirichlet boundary condition u = on dD instead of 
the Neumann boundary condition du/du = on dD, one can drop the assumption on uj 
and the Gaussian curvature of dD. See [10] for this result. Thus we propose 
Open problem 2.1. Remove the curvature condition in Theorem 2.1. 

A sketch of the proof of Theorem 2.1. Let u solve (1.1) with / = v\$n and set w = 
u — v inf2 \ D. The w satisfies (1.2). We have three lemmas. 

Lemma 2.1. There exists a positive constat C{k) such that for all uj G S 2 , r > 

2r 2 f e 2rx-u dx _ k 2 j , | 2(fa < f , A _ Ac )( v | an ) . ^ < C(A;)(r 2 + A; 2 ) / e 2 ™'"^. 
Jd Jq\d Jan ' jd 

This is a consequence of the representation formula (1.3) and the estimate HHIif^nVD) — 
C{k)\\v\\ H i( D) . 



Lemma 2.2. 



liminfe' 2r ^ {w) r 2 / e 2rx ^dx > 0. 

T >00 J D 



The proof of this lemma can be done by slicing D with the planes x ■ uj = ho{oo) — s 
with < s « 1. 

Lemma 2.3. Assume that the set {x G dD \ x • uj — ho(oj)} consists of the only one point 
and the Gaussian curvature of dD doesn't vanish at the point. Then 

_ \w\ 2 dx 
lim - = 0. 

r ^°° 2r 2 / e 2rx ^dx 

JD 

From Lemmas 2.1, 2.2 and 2.3 one knows that there exist positive constants C\, C 2 
and t > such that for all r > r 

All the statements in Theorem 2.1 now follows from these estimates. 

Finally we describe the outline of the proof of Lemma 2.3. One can find p G H 2 (Q\D) 
such that (A + k 2 )p = w in Q \ D, p = on <9f2 and dpjdv = on dD. From the Sobolev 
imbedding and the estimate Ibllfp/m™ < C(k) ||HI.L 2 (nYD) we have: \p(x) — p(y)\ < 
C(k)\x -yl^WwW^^ and sup \p(x)\ < C^HHL^vd)- 

Let x be the point in the set {x G <9-D I x ■ uj — h D (uj)}. Since / (dv/dv)dS(x) = 

JdD 

—k 2 / vdx, one can write 

JD 

( _\w\ 2 dx = — [ p^-dS(x) = f \p(x ) — p(x)}^-dS(x) + k 2 p(x ) [ vdx. 
Jn\D JdD dv JdD dv Jd 

From these one gets 

J \w\ 2 dx < C{k) (V2t 2 + k 2 \x - x\^ 2 e TX - u dS(x) + e TX "dx^j IMIl^vd) 



11 



and this thus yields 

\w\ 2 dx < C{k) | (tJ 9d \x - x\ 1/2 e TX - w dS(x)^ + (^J D e TX0J dx 
The Schwarz inequality yields 

(^j D e TXW dxy < \D\ J e 2rx - u dx. 
Thus from this and Lemma 2.2 one knows that it suffices to prove 

lim re - ThD{uj) [ \x - x\ 1/2 e TX - u dS(x) = 0. 

T — >°° JdD 

In fact, one gets 

re- ThD ^ [ b - x^e^dSix) = O^ 1 / 4 ). 

JdD 

This is proved by using a localization at x and a local coordinates at the point. 



2.2 The enclosure method with a single incident plane wave 

The idea started with considering an inverse boundary value problem for the Laplace 
equation in two dimensions in [8]. Five years later in [15] the idea was applied to an 
inverse obstacle scattering problem in two dimensions. The problem is to reconstruct a 
two dimensional obstacle from the Cauchy data on a circle surrounding the obstacle of 
the total wave field for a single incident plane wave with a fixed wave number. 

In this subsection we assume that D is polygonal, that is, D takes the form DiU- • -UD m 
with 1 < m < oo where each Dj is open and a polygon; Dj D Dy = if j ^ f. 

The total wave field u outside obstacle D satisfies 

011 

Au + k 2 u = ® inR 2 \D, — = on<9£> 

ov 

and the scattered wave w = u — e lkx ' d with k > and d E S 1 satisfies the outgoing 
Sommerefeld radiation condition lim Jridw Idr — ikw) = 0, where r = \x\. 

r — >oo 

Let B R be an open disc with radius R satisfying D C Br. We assume that Br 
is known. Our data are u and dujdv on 8Br. Let u and u 1 - be two unit vectors 
perpendicular to each other. Set z = rco + iy/r 2 + k 2 ^ 1 - with r > and v(x; z) = e x ' z . 
Recall hn{uj) = sup x£D x ■ u. 

Theorem 2.2. Assume that the set dD n {x e R 2 | x ■ u = hr>(uj)} consists of only one 
point. Then the formula 



lim — log 

r >oo j- 



f ( Ou Ov \ 

L R w {x]z) -Tv {x]z)u ) dS{x) 



h D (u), 



is valid. Moreover, we have: 
ift> h D {uj), then 



lim 

T >00 



L R \^ e ~ Ttv ^ z) - e ~ T V^ z)u ) dS{x) 



0; 
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if t < hoioj), then 



lim 

r >oo 



LM e ~ T,viX]Z) - e ^ {x '- z)u ) dS{x) 



OO. 



Sketch of the proof. The one of key points is: introducing a new parameter s instead 
of t by the equation s = \Jt 2 + k 2 + r, we obtain, as s — > oo the complete asymptotic 
expansion 



L 



^(x; z) - ^(x; z)u dS(x)e- WT ' +kix °- u - ThD M ( 2 .1) 

Here the A n describes the singularity of tt at a corner and in this case explicitly given by 
the formula X n — (n — l)n/Q, where denotes the outside angle of D at x G <9-D n {x e 
R 2 | x • w = /i.d(u;)} an d thus satisfies 7r < < 27r; fT n are constants depending on X n , 
u> and shape of D around x ; a 2 , «3, • • • are the coefficients of the convergent series 
expansion of u with polar coordinates at a corner: 

oo 

u(r,6) = aiJo(kr) + Oi n J\ n (kr) cos\ n 6, < r << 1, < 9 < 0. 

n=2 

Now all the statements in Theorem 2.2 follow from (2.1) and another key point: 3n > 2 
a n K n 7^ 0. This is due to a contradiction argument. Assume that the assertion is not 
true, that is, Vn > 2 a n K n = 0. 

First we consider the case when Q/ir is irrational. In this case we see that Vn > 
2 fT n 7^ 0. Thus a n = and this yields u(r,9) = a\Jo(kr) near a corner. Since this 
right-hand side is an entire solution of the Helmholtz equation, the unique continuation 
property of the solution of the Helmholtz equation yields u(x) = aiJo(k\x — xo\) in H 2 \D. 
However, we see that the asymptotic behaviour of this right-hand and left-hand sides are 
completely different. Contradiction. 

Next consider the case when Q/n is a rational. By carefully checking the constant 
K n we know that for each n > 2 with K n = the X n becomes an integer. From the 
assumption of the contradiction argument one knows if n satisfies K n ^ 0, then C n = 0. 
Thus we have the expansion 

u(r, 9) =J2 Cn 3 J\ nj (kr) cos X nj 9, 

rij 

where rij > 2 satisfy K n = 0. Since X n . is an integer and A n = (rij — l)ir, from this 
right-hand side one gets: for all r with < r << 1 du/d9(r, it) = du/d9(r, — n) = 0. 
Then a reflection argument ([1]) yields that this is true for all r > 0. However, from this 
together with the asymptotic behaviour of Vtt one can conclude that incident direction d 
has to be parallel to two linearly independent vectors which are directed along the lines 
9 = 7i and 9 = — tt. Contradiction. 
Remarks are in order. 

• In Theorem 2.2 one uses the Cauchy data on the circle surrounding the obstacle as the 
observation data. However, du/du on Br can be calculated from u on 8Br by solving an 
exterior Dirichlet problem for the Helmholz equation. 
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• In [18] a similar formula has been established by using the far field pattern Fd(<p, d; k), 
(p E S 1 of scattered wave w = u — e lkx ' d for fixed d and k which determines the leading 
term of the asymptotic expansion of w at infinity in the following sense 



Moreover, therein instead of volumetric obstacle, similar formulae for thin sound-hard 
obstacle (or screen ) also have been established with two incident plane waves. 

• In [37] the numerical testing of a method based on results in [14, 15, 18] has been 
reported. 

• It would be interesting to consider the case when the total wave u satisfies the equation 
V • 7VW + k 2 u = in R 2 where 7(2;) = 1 for x G R 2 \ D and 7(2;) = Aj for x e Dj, 
j = l,---,m; each Aj are positive constants and Aj ^ 1. The author thinks that this 
case becomes extremely difficult because of the complicated behaviour of u at a corner. 
However we propose 

Open problem 2.2. Establish Theorem 2.2 for u above. 

See [11] for k = and [17] for the equation V • 7VM + k 2 = 0. 

• For recent applications of the enclosure method with a single measurement for a system 
arising in linear theory of elasticity we have [24, 25, 26]. However, their extension to the 
elastic wave with a single incident plane wave remains open. It is a challenging problem 
to be solved. 

3 Inverse obstacle scattering problems with dynam- 
ical data over a finite time interval 

Previously we considered only the stationary or time harmonic problem. In this section 
we consider how one can use the data over a finite time interval to extract information 
about the location and shape of unknown obstacles. In [28, 39, 40] some uniqueness 
results have been established, however, it seems that mathematically rigorous study of 
the reconstruction issue in this type of problem has not been paid much attention. Note 
that: there are some results [31, 32, 33] in the context of the Lax-Phillips scattering 
theory, which give the convex hull of an unknown obstacle, however, the data are taken 
from t = to t = 00. 

The purpose of this section is to introduce a new and simple method in [23] which is 
an application of the idea developed in [22, 27] and employs the data over a finite time 
interval on a known surface surrounding unknown obstacles. 

3.1 New development of the enclosure method 

In order to explain the basic idea, in this subsection we present an application to the 
one-space dimensional wave equation which is taken from Appendix B in [22]. Let a > 
and c > 0. Let u — u(x, t) be a solution of the problem: 




2 U tt = Uxx in]0, a[x]0, T[, cu x (a,t) = fort e]0, T[, 



u(x,0) = 0, u t (x,0) = in]0, a[. 
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The quantity c denotes the propagation speed of the signal governed by the equation. 
Inverse Problem 3.1. Assume that a is unknown. Extract a from u(0,t) and u x (0,t) 
for < t < T. 

Theorem 3.1. Let u x (0,t) e L 2 (0, T) satisfy the condition: there exists a real number fi 
such that 

rT 



lim inf r M 



/ u x (0,t)t 
Jo 



- Tt dt 



> 0. 



Let T > 2a/ c and v(x, t) = v(x, t;r) — e r( - x / c+t \ Then the formula 



lim — log 

T >00 f 



[ T (-cv x (0, t)u(0, t) + cuJO, t)v(0, t)) dt 
Jo 



-2a/ c, 



(3.1) 



(3.2) 



is valid. 



Some remarks are in order. 

• The v satisfies the wave equation (l/c) 2 v tt = v xx and satisfies: if x + ct > 0, then 
v (x, i) — > as t — > oo; if x + ct < 0, then v(x, t) — > +oo as r — > oo. 

• The quantity 2a/ c coincides with the travel time of a signal governed by the wave 
equation with propagation speed c which starts at the boundary x = and initial time 
t = 0, reflects another boundary x = a and returns to i = 0. Thus the restriction 
T > 2a /c is quite reasonable and does not against the well known fact: the wave equation 
has the finite propagation property. 

• The condition (3.1) ensures that u x (0,t) can not be identically zero in an interval 
]0, T"[c]0, T[. Therefore surely a signal occurs at the initial time. However, it should be 
emphasized that the formula (3.2) makes use of the averaged value of the measured data 
with an exponential weight over the observation time. This is a completely different idea 
from the well known approach in nondestructive evaluation by sound wave: monitoring 
of the first arrival time of the echo, one knows the travel time. 

A sketch of the proof of Theorem 3.1. Introduce the function w by the formula 

rT 

w(x) — w(x;t) — / u(x,t)e Tt dt, < x < a. 
Jo 



It holds that 



2 

C W 



t 2 w = e TT (u t (x,T) + tu(x,T)) in]0, a[, cw'(a) = 0. 



Then, this together with integration by parts gives the expression 



2ar/c 



( (-cvJO, t)u(0, t) + cu x (0, t)v(0, t)) dt 
Jo 



Tw(a)e aT/c - c^e-^-t 2 



a ^ (\u t {i,T) + ru{i,T))e-^/ c di. 
Jo 



Now (3.2) can be checked by studying the asymptotic behaviour of this right-hand side 
with the help of the expression 



w(a) 



2cw'(0) 



-tT 



t e 



ar/c 



-ar/c 



) 4 



ar/c 



-ar/c 



) 



X 



{ j o a M£,T) + ru{i,T))e-^' c di + J\u t {i,T) + Tu{t,T))e^ c dt 
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together with (3.1). 

The proof presented here heavily relies on the spaciality of one-space dimension. In 
[27] we found another method for the proof which works also for higher space dimensions 
and applied it to a similar problem for the heat equation. In the following two subsections 
we present further applications of the method to the wave equations. 

3.2 Sound-hard obstacle 

Let D C R 3 be a bounded open set with smooth boundary such that R 3 \ D is connected. 
Denote by v the unit outward normal to 3D. Let < T < oo. 

Given / G L 2 (R 3 ) with compact support satisfying supp / fl D — let u — u(x,t) 
satisfy the initial boundary value problem: 

f)n 

d 2 t u - Au = in(R 3 \ D)x ]0, T[, — = ondDx ]0, T[, 

u(x, 0) = 0, d t u(x, 0) = f(x) in R 3 \ D. 

Let Vl be a bounded domain with smooth boundary such that D C VL and R 3 \ Vt is 
connected. Denote by the same symbol v the unit outward normal to dQ. 

The dQ is considered as the location of the receivers of the acoustic wave produced by 
an emitter located at the support of /. In this section we consider the following problem. 
Inverse Problem 3.2. Assume that D is unknown. Extract information about the 
location and shape of D from u on <9f2x]0, T[ for some fixed known f satisfying supp/n 
Q = and T < oo. 

Note that u in (R 3 \ f2)x]0, T[ can be computed from u on <9f2x ]0, T[ by the formula 

u = z in(R 3 \n)x]0, T[ (3.3) 
where z solves the initial boundary value problem in R 3 \ f2: 

d 2 t z - Az = in(R 3 \n)x]0, T[, z = u on9fix]0, T[, 

z(x, 0) = 0, d t z(x, 0) = f(x) in R 3 \ U. 

Thus the problem can be reformulated as 
Inverse Problem 3.2'. Extract information about the location and shape of D from u 
in (R 3 \ Q) x]0, T[ for some known / satisfying supp / fl Q = and T < oo. 

Now we state the result. Let B be an open ball with B fl Q = 0. Choose the initial 
data / G L 2 (R 3 ) in such a way that: 

(11) f( x ) = a.e. x G R 3 \B; 

(12) there exists a positive constant C such that f(x)>Ca.e. x G B or — f(x) > C a.e. 
x G B. 

Set 

rT 

w(x;t)= / e~ Tt u(x,t)dt, xeK 3 \D, r > 0. 
Jo 

Our result is the following extraction formula from w and dwjdv on dQx ]0 T[ which can 
be computed from the data u in (R 3 \ Q)x ]0, T[. 



16 



Theorem 3.2. Let r > and v G /J 1 (R 3 ) be the weak solution of 

(A - t 2 )v + f(x) = inR 3 . (3.4) 
If the observation time T satisfies 

T > 2dist(D, B) — dist(Q, B), (3.5) 
then there exists a r > such that, for all r > r 



and the formula 



Jan \dv dv ) 

1 f I dv dw \ 
lim — log/ (3.6) 

r _^^2r B Jan\dv dv J v ' ; ' v ; 



Some remarks are in order. 
• The v is unique and is given by the explicit form 

1 r e~ T)[X ~ y]{ 



v(x;t) = — I 7 r/(j/)dy, a: G R 3 . 

47r 7s |x — y\ 



• The quantity dist (D, B) + y \dB\/Aix coincides with the distance from the center of B 
to D and thus (3.6) yields the information about d£>(p) for a given point p in R 3 \ Q. 

• It is easy to see that 2dist (D, B) - dist (ft, 5) > Z(9S, dD, dtt), where Z(9S, 9D, dfl) = 
inf — y| + \y — z\ \ x G dB ,y G z G <9fi}. This is the minimum length of the broken 
paths that start at x G dB and reflect at y G <9-D and return to z G <9f2. Therefore (3.5) 
ensures that T is greater than the first arrival time of a signal with the unit propagation 
speed that starts at a point on dB at t = 0, reflects at a point on dD and goes to a point 
on dQ. 

The main part of the proof of Theorem 3.2 is to show that 

lim inf T ^ dist { d,b) f (d^ w _ dw\ dS > Q (3 J) 

T — >°° Jan \dv dv ) 

It is a consequence of the following representation formula which corresponds to (1.3) and 
the estimate for v: 

/ 

Jan \dv dv J 

= ( \Vv\ 2 dx + T 2 ( \v\ 2 dx+ [ \V(w -v)\ 2 dx + r 2 [ \w - v\ 2 dx 
Jd Jd Jr3\d Jr 3 \d 

+e~ rT ( _(w -v)(d t u(x,T) + ru{x,T))dx - e~ rT [ _(d t u(x,T) + ru{x,T))vdx; 
Jr 3 \d Jn\D 

liminfr 6 e 2Tdist ( D < B ) / \v\ 2 dx > 0. (3.8) 
r— ^oo Jd 

Note that the precise values of 4 and 6 of r 4 in (3.7) and r 6 in (3.8), respectively are not 
essential. 
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3.3 Penetrable obstacle 



The method in the former subsection can be applied to a more general case. Given 
/ G L 2 (R 3 ) with compact support let u = u(x,t) satisfy the initial value problem: 



where 7 = 7(2;) = (7^(2?)) satisfies: for each i,j = 1,2,3 7y(x) = 7^(2;) G L°°(R 3 ); there 
exists a positive constant C such that "f(x)£ ■ £ > C|£| 2 for all £ G R 3 and a. e. x G R 3 . 

We assume: there exists a bounded open set D with a smooth boundary such that 
7(2;) a.e. x G R 3 \ D coincides with the 3x3 identity matrix J 3 . Write h(x) = 7(2;) — J 3 
a.e. x G -D. 

Our second inverse problem is the following. 
Inverse Problem 3.3. Assume that both D and h are unknown and that one of the 
following two conditions is satisfied: 

(Al) there exists a positive constant C such that —h(x)£ ■ £ > |£| 2 for all £ G R 3 and a.e. 
x G D; 

(A2) there exists a positive constant C such that h(x)£ ■ £ > |£| 2 for all £ G R 3 and a.e. 
x e D. 

Let be a bounded domain with smooth boundary such that D C Q. Extract information 
about the location and shape of L> from u on 9Qx]0, T[ for some fixed known f satisfying 
supp / l~l H = and T < 00. 

Note that u in (R 3 \ Q) x]0, T[ can be computed from u on 5i7x]0, T[ by the exactly 
same formula as (3.3) and thus the problem can be reformulated again as 
Inverse Problem 3.3'. Extract information about the location and shape of D from u 
in (R 3 \ Q)x ]0, T[ for some known / satisfying supp / fl f2 = and T < 00. 

Now we state our second result. 

Theorem 3.3. Assume that 7 satisfies (Al) or (A2). Let f satisfy (II) and (12) in 
subsection 3.2 and v be the weak solution of (3.\). LetT satisfies (3.5) and w be given by 



with solution u of (3.9). If (Al) is satisfied, then there exists a r > such that, for all 
r>r 



d 2 t u - V • 7VW = inR 3 x ]0, T[, 



(3.9) 



u(x,0) = 0, d t u(x,0) = f(x) inR 3 , 





if (A2) is satisfied, then there exists a r > such that, for all r > r 




In both cases we have 




dist(D,B). 
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The key points for the proof are an estimate for Vi> similar to (3.8) and the following 
two representation formula: 



/ {(Vi> • v)w — (7VW • v)v\dS = — I hVv ■ Vvdx 
Jan Jd 

+ / 7V (u> — v ) ■ V(w — v )dx + r 2 / \w — v\ 2 dx 
+e~ rT / (d t u(x, T) + tu(x, T))(w — v)dx — e~ rT / (d t u(x, T) + ru(x, T))vdx; 

— / {CVv ■ u)w — ('yVw ■ u)v}dS — / ■ Vtudi 

Jan Jd 

+ / V(i> — w) ■ V(f — w)dx + t 2 / |t> — w\ 2 dx 

JR3 JR3 

tu(x, T))vdx. 

J R,^ J Q 



4 Summary and further research direction 

In this paper we presented: past applications of the probe and enclosure methods to 
inverse obstacle scattering problems with a fixed wave number and related open problems; 
recent applications of the enclosure method to inverse obstacle scattering problems with 
dynamical data over a finite time interval. 

In particular, in Section 3 we presented a new and simple method in [23] for a typical 
class of inverse obstacle scattering problems that employs the values of the wave field over 
a finite time interval on a known surface surrounding unknown obstacles as the observation 
data. The wave field is generated by an initial data localized outside the surface and its 
form is not specified except for the condition on the support. The method explicitly yields 
information about the location and shape of the obstacles more than the convex hull. 

It would be interesting to apply the method presented in Section 3 to other time 
dependent problems in electromagnetism(e.g., subsurface radar [4], microwave tomography 
[41]), linear elasticity, classical fluids etc.. Those applications belong to our future plan. 
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